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Summary 

COSP is a publicly available computer program that models 
the cyclic oxidation weight gain and spallation process. Inputs 
to the model include the selection of an oxidation growth law 
and a spalling geometry, plus oxide phase, growth rate, spall 
constant, and cycle duration parameters. Output includes weight 
change, the amounts of retained and spalled oxide, the total 
oxygen and metal consumed, and the terminal rates of weight 
loss and metal consumption. The present version is Windows 
based and can accordingly be operated conveniently while 
other applications remain open for importing experimental 
weight change data, storing model output data, or plotting model 
curves. Point-and-click operating features include multiple 
drop-down menus for input parameters, data importing, and 
quick, on-screen plots showing one selection of the six output 
parameters for up to 10 models. A run summary text lists various 
characteristic parameters that are helpful in describing cyclic 
behavior, such as the maximum weight change, the number of 
cycles to reach the maximum weight gain or zero weight change, 
the ratio of these, and the final rate of weight loss. The program 
includes save and print options as well as a help file. Families 
of model curves readily show the sensitivity to various input 
parameters. The cyclic behaviors of nickel aluminide (NiAl) 
and a complex superalloy are shown to be properly fitted by 
model curves. However, caution is always advised regarding 
the uniqueness claimed for any specific set of input parameters. 

Introduction: Model Overview 

COSP for Windows (WinCOSP is the abbreviated program 
file name) is an adaptation of original cyclic oxidation models 
developed at NASA Glenn Research Center over a period of 
25 years (refs. 1 to 8). Some papers documenting many of the 
precepts and findings are listed in these references and related 
bibliographies (refs. 9 to 15). Original mainframe versions were 
written in Fortran, and later a version for personal computers 
was written in Modula-2 and described in a complete archival 
treatise (ref. 7). It is hoped that this Windows version will further 
increase its user-friendly attributes and broaden its base of users. 


Basically, the program calculates a cyclic oxidation weight 
change curve by using some elementary inputs and assump- 
tions regarding model scale growth and spalling behavior. 
Along with the sample weight change, other factors germane to 
the cyclic process — oxygen consumed, metal consumed, scale 
weight, spall weight — are all part of the calculations. Key 
descriptive parameters, such as number of cycles to maximum 
weight gain, cycles to cross zero weight change, final weight loss 
slope, and the plateau in oxygen gain, are all highlighted as part 
of the text output. 

The mechanics of the calculation are briefly described by 
way of introduction. A type of oxide is chosen (e.g., alumina, 
AI2O3) in order to determine a stoichiometric constant S f , or the 
weight ratio of oxide to oxygen in the scale. This constant 
relates the amount of oxide weight (metal + oxygen) to that of 
oxygen in the scale (e.g., 2.124 for AFO3). A growth model 
may be chosen, but typically parabolic rates are assumed, where 
the specific weight change A WfA is given by equation (1). 
Options for power law (eq. (2)) and logarithmic grow th (eq. (3)) 
are also included: 


■4*; 

II 

rj 

<1 

(1) 

(A W/A) m = kt 

(2) 

AW/A = ln b/ + c) 1/m ] 

(3) 


where k p is the parabolic rate constant, k is the grow th rate for 
the other models, t is time, c is a secondary growth constant, and 
m is a growth exponent. 

The growth model is used to calculate the amount of scale 
grown at temperature during one cycle. The amount of spalla- 
tion that occurs during the cooling portion of the cycle is 
governed by the spall constant Qq. The fraction spalled F s of the 
existing scale weight W' is given by 

F s ={Qo){w;) a (4) 


NASA/TP— 2002-21 1 108 


1 


where a i$ the spall exponent, a constant usually assumed to be 
1.0 (as found in some independent tests of spall amount versus 
scale weight, refs. 2, 4, and 7). It is important to note that Q 0 is 
not strictly a spall fraction or percentage. The remaining amount 
of oxide after C cycles W r c and the effective time t cff c that 
corresponds to the isothermal growth of this amount of 
oxide are then calculated. The scale present during the next 
heating cycle W' c+ j can then be determined on the basis of 
the growth model, but now by using the time + At, where At 
is simply the cycle duration. This sequence of calculations is 
seen by inspection of equations (5) to (7) below. Repetition of 
these calculations for a preset number of cycles generates the 
scale weight and allows calculation of the entire cyclic oxida- 
tion curve. 

w r = \v;-F s w; (5) 

Kc = S r (^Af) 1/2 (6) 

w;, c+ i=^[^(/ effiC +A/)] 1/2 (7) 

In addition to Q 0 as an input for spalling, the spallation 
“geometry” must also be specified. The spallation geometry 
describes whether the material is removed as a uniform layer off 
the top of the entire scale or whether it is removed in discrete 
areas as “chunks” of specified thickness. For many instances 
the uniform spalling layer provides a reasonable approximation 
of the average behavior of systems that do not exhibit primarily 
interfacial spallation. Nonuniform or discrete spallation areas 
are more precisely handled by the Monte Carlo bimodal calcu- 
lations. Within this option the sample surface is divided into a 
specified number of equivalent segments. The total spall frac- 
tion F s is given by F = P{F spa nk where P is the probability 
and tf spa |j is the depth fraction (specified from the input menu) 
of the spalled segment for a bimodal distribution of spall 
geometry. A random number generator (0 to 1 .0) is used to 
provide a random value RN for each segment, and when 
(RN) < P, that segment is tagged for spallation of the fractional 
depth P spa ]] on that cycle. The total spall fraction of the scale is 
summed over all the segments and is equal to F av(7 = Q 0 ( W' ) a . 
The specific scale measures for each segment are then cata- 
loged and updated each cycle for the entire oxidation curve. 

One final option is to use a rectangular distribution of spall 
probability (ref. 7), but in practice it is found to be indistin- 
guishable from the uniform spalling case. The Windows pro- 
gram also allows for importing tables of experimental specific 
weight change versus time data, so that model curves can be 
directly compared with experimental results. 

This report reviews the overall findings of the original COSP 
model as presented previously in references 3 to 8, with 
reference 7 offering the most complete development of the 
theory and critical analysis of the model results. These findings 


are presented in the format of what the COSP for Windows user 
experiences in stepping through the main windows and model 
options. Finally, having established that background, a few 
examples are given that show how the models can be used to fit 
the behavior of well-known, high-temperature, oxidation- 
resistant materials. 


COSP for Windows requires Windows95 or later versions. 
Three windows are always accessible by large point-and-click 
buttons: the Plot window, the Results window, and the 
Experimental data window, whose major features are de- 
scribed below. (The specific window titles, menu lists, and 
buttons are displayed here in boldface type.) 


Under the Plot window, the Setup dialog window is used to 
define the model parameters and plot the curv es. It appears when 
the program is started, using the initialization default param- 
eters (Af0 3 , k = 0.01 mg-/cm 4 hr; uniform spall with 
Qq = 0.01, a = 1, and At = 1 hr). The basic inputs are listed 
below: 

1 . Oxide (Stoichiometric constant), S c 

2. Growth law: parabolic, power law, or logarithmic 

a. Parabolic: (Specific weight gain) 2 = k p t 

b. Power law: (Specific weight gain) m = kt 

c. Logarithmic: Specific weight gain = ln[(A? + c) 1/,w ] 
Rate constants, k p (mg 2 /cm 4 hr) or k 

Growth exponent, m 
Secondary growth constant, c 

3. Spall constant, Q 0 

4. Spall exponent, a 

5. Spall case: 

a. Uniform: F s = (f> 0 )( W' ) a 

b. Monte Carlo: Bimodal F assigned from 

F , Vg = m W r) a =P(R spM l 

Segment spalls if P* < P, where P* is assigned by a 
random number generator for each segment. 

6. Number of segments 

7. Number of cycles, C 

8. Time per cycle, At , hr 

An instant plot is shown in the graphic display field when the 
Okay button is clicked. The Add another model button allows 
the user to simultaneously display the results for up to 10 
different models. Different model outputs for the y axis may be 
selected (e.g. , Cumulative oxygen uptake, Cumulative weight 
of metal consumed, Fraction of oxide spalled, Weight of 
retained oxide, etc.) or, for the x axis, Time (in hours) or 


Plot Window 


Basic Operating Features 
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Cycles. The resulting plot may be printed directly to any 
Windows-compatible printer. A short input listing will be printed 
for each existing model. The model parameters may be saved 
directly in a COSP for Windows folder for easy retrieval and 
future use in the program. 

Results Window 

A tabular display of all the calculated output data for the last 
model selected is presented in the Results window, as is a text 
summary of the input parameters and a description of the curve 
(see output list below). They may both be copied and pasted into 
spreadsheet or text applications. The model menu allows one 
to select and display the data for other models or to delete all the 
current models. The outputs are as follows: 

1. Summary (text): 

a. Cycles to maximum weight change, C m 

b. Maximum weight change, (AW/A) max 

c. Cycles to zero weight change, C 0 

d. Ratio C^/Cq, or inverse 

e. Final slope (negative weight loss rate) 

f. Final metal consumption rate 

g. Final effective time, t eff 

2. Spreadsheet (table): 

a. Cumulative weights for 

Oxygen uptake 
Metal consumed 
Spalled oxide 

b. Weight of retained scale 

c. Fraction of oxide spalled, F s 

d. Net weight change, AW/ A 

e. Effective time, r eff 


Experimental Data Window 

Sets of experimental weight change data may also be imported 
by using a Paste from clipboard command in the Experimen- 
tal data window. These data will be then plotted in the Plot 
window. This feature can be used to directly compare experi- 
mental oxidation curves with the results from the calculated 
models. The total number of experimental and model curves 
cannot exceed 10. 


General Model Trends 

Baseline Behavior 

In this section some basic trends in the cyclic oxidation curve 
are reexamined as a function of all the input parameters. A 
standard baseline for comparison is offered assuming an A^O^ 
scale with a parabolic growth rate of 0.01 mg 2 /cm 4 hr (represen- 
tative of 1200 °C growth), uniform layer spallation with a Q 0 
of 0.002 cm 2 /mg, a cycle duration of 1 hr, and a test length 
of 1000 cycles. Many of the outputs for this model are plotted 
against time (cycles) in figure 1. The cyclic sample weight 
change curve reaches a maximum weight gain of about 
0.66 mg/cm 2 at 85 cycles and reaches a negative weight change 
at 279.2 hr with a ratio C ;w /Cq of 0.30. The final slope in 
AW/A (linear weight loss rate) is 0.0053 mg/cnrhr. As this 
linear rate is approached, so too are the linear rates of oxygen 
(I0 2 ) and aluminum (LAY) consumption as well as that of total 
oxide spalled (IlV spa11 ). (It is useful to note that the final slope 
of the weight change curve is equivalent to the final rate of 
metal consumed (ref. 7).) Figure 1(a) shows that the net 
weight change curve is equal to the total oxygen reacted less the 



Number of 1 -hr cycles, C 



Figure 1 .—Typical outputs from standard baseline COSP cyclic oxidation model. Parabolic growth and uniform spalling of 
Al 2 0 3 scale with parabolic rate constant k p , 0.01 mg 2 /cm 4 hr; spall constant Q 0 , 0.002 cm 2 /mg; and cycle duration At, 

1 hr. (a) Net sample weight change AW/A related to total amount of oxygen (SO 2 ) and aluminum (SAI) consumed and 
amount of oxide spalled (2 W spa j|). (b) Net sample weight change related to existing scale weight W r and scale fraction 
spalled F s each cycle. 
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amount of oxide spalled (X0 2 - Z^ spa]I ). Figure 1 (b) shows that 
the retained scale and consequently the fraction of scale spalled 
in each cycle approach a constant level at the same time that 
linear rates of weight loss, material consumed, and total oxide 
spalled are approached. These parameters achieve better than 
92 percent of the final limits when the number of cycles C = C Q 
and 99.9 percent of same when C = 2 C 0 . This steady-state 
behavior suggests that some extrapolation from relatively short- 
time testing is warranted from the model, provided that no 
mechanism changes take place. 

Parametric Trends 

The sets of weight change curves displayed in figures 2 to 6 
represent major trends with the various input parameters. A 
baseline weight change curve (from fig. 1) is offered for com- 
parison in each figure. 

Stoichiometric constant , — In practice any variation in scale 
phase usually entails a consequential change in growth rate at 
any given temperature. However, for the purpose of showing the 
singular effect of varying the stoichiometry and composition of 
the scale alone, figure 2 displays a family of cyclic oxidation 
curves for alumina (A1,0 3 )* chromia (Cr 2 0 3 ), nickel aluminate 
spinel (NiAI 2 0 4 ), nickel oxide (NiO), and nickel tantalate 
(NiTa^O^). It is apparent that, in order of increasing S c , the 
cycles to reach maximum weight gain C and the cycles 
to cross zero Cq steadily decrease, whereas the terminal 
(negative) weight loss slope steadily increases. This result is 
the effect of altering only the scale (and spallation layers) to 
those with heavier metal (weight) contents with respect to the 



Figure 2. — Effect of oxide type (stoichiometric constant) 
on baseline COSP model curves for common (Ni-Cr-Af) 
oxides. Parabolic growth and uniform spading case with 
parabolic rate constant k p , 0.01 mg 2 /cm 4 hr, spall 
constant Qq, 0.002 cm 2 /mg; and cycle duration A f, 1 hr. 



Number of 1 -hr cycles, C 




Figure 3. — Effect of parabolic rate constant, spall constant, 
and cycle duration on standard baseline COSP model 
curves for parabolic growth and uniform spalling of AI2O3 
scale. Parabolic rate constant k pt 0.01 mg 2 /cm 4 hn spall 
constant Qq, 0.002 cm 2 /mg; and cycle duration At = 1 hr. 
(a) k p varied, (b) Qq varied, (c) At varied. 
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oxygen (weight) content. (The use of identical growth rates is 
unrealistic and is done here only for the sake of illustration.) 
Simple combined oxide stoichiometries can be used to repre- 
sent multiple phases in the scale, such as NiA^O^ to represent 
NiO + A1 2 0 3 . If a complex stoichiometry is needed and cannot 
be easily represented by a chemical formula, it is possible to 
enter any stoichiometric constant. For example, the arithmetic 
average of two constants may be entered to represent a scale 
consisting of a 50-50 percent mixture (by weight) of the two 
oxides. 

Rate constant . — The effect of varying the rate constant k : in 
the case of parabolic growth and uniform layer spallation, 
keeping Q 0 at 0.002 and cycle duration at 1 hr, is shown in 
figure 3(a) for an AI 2 0 3 scale. Increasing k p increases the 
maximum weight gain and the rate of final weight loss and 
decreases the cycles to both maximum and zero weight gain. 
The ratio CJC Q remains essentially constant at 0.305 
to 0.310. 

Spall constant . — Figure 3(b) shows the effect of increasing 
the spall constant Q 0 for the case of parabolic growth and 
uniform layer spallation, fixing k p at 0.01 mg 2 /cm 4 hr. Also 
shown for comparison are the two extremes of (a) zero spalling 
for Q 0 = 0.0 and (b) total spallation. For the case of total 
spallation in each cycle 

F s = 1.0 (8) 

From equation (4) one obtains 

^=( 0 o )(w;) a = i.o (9) 

and from equation (6) 


/ \ 1 / 2 

w; = s c (k p At) 

one derives that 

C?0.max = 7 * ,|/ 2 - = 4 ' 7074 < 10 > 

ScMO 

for this case of A1 2 0 3 . Between the two extremes, as Q 0 
increases, the maximum weight gain steadily decreases, as does 
the number of cycles to maximum and zero weight gain. The 
ratio C^/Cq remains at about 0.31, where applicable. The final 
rate of weight loss increases with Qq. 

Cycle duration . — Finally, figure 3(c) shows the effect of 
cycle frequency (or duration). Here it is important to note 
whether the outputs are plotted versus cycle number or oxida- 
tion time. In this plot the maximum weight gain decreases with 
increasing cycle frequency, as do the number of hours to 
maximum and zero weight change. The ratio CJC 0 remains 
constant between 0.304 and 0,308. The final rate of weight loss 
increases with cycle frequency. 

Growth law. — In most cases using a parabolic growth law is 
adequate for cyclic oxidation modeling. Deviations from para- 
bolic kinetics are generally minor or require long oxidation 
times to become evident. Also, essentially parabolic curves can 
be fitted very nicely by adjusting both parameters ( k and m) in 
the power growth law or logarithmic growth law models. This 
section illustrates some trends with growth parameters for these 
two other models. The baseline model is for A1 2 0 3 with a spall 
constant Q 0 of 0.002 cm 2 /mg and a cycle duration A/ of 1 .0 hr. 

Figure 4 shows the families of power growth law curves for 
varying k and m. In figure 4(a) increasing k for the cubic 




0 200 400 600 800 1000 

Number of 1-hr cycles, C 


Figure 4.— Effect of power law growth kinetics on baseline COSP model curves for uniform spalling of AI 2 C>3 scale. Spall 
constant Q 0 , 0.002 cm2/mg; cycle duration Af, 1 hr. (a) Growth rate k varied for growth exponent m of 3. (b) m varied for k 
of 0.01. 
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(m = 3) case produces effects analogous to those shown in 
figure 3(a) for increasing k in parabolic growth. Again the ratio 
C n /C 0 remains constant but now lies between 0.223 and 0.226. 

In figure 4(b) subparabolic behavior (m > 2) lessens the final 
weight loss rate and shortens the time to reach maximum 
weight without appreciably shortening the time to reach zero 
weight change. Such models may be useful where fast-growing 
initial transient oxides give way to slower growing more stable 
scales or where oxide grain growth decreases the relative area 
of short-circuit pathways with time. Conversely, linear behav- 
ior (m = 1 .0) lengthens the time to reach maximum weight and 
produces greater weight loss rates, as might be expected for a 
faster growing scale. The ratio C m /C 0 is 0.47 for linear (m = 1 ), 
0.30 for parabolic (m = 2), and 0.22 for cubic (m = 3) behavior. 

Figure 5 shows families of logarithmic growth law curves for 
varying the rate constant k and the growth exponent m. In fig- 
ure 5(a) the secondary growth constant c, which is the argument 
for the logarithm at t = 0, is set to 1 .0 in order to yield zero- time 
weight change of 0,0. And m is set at 1.0 so that the weight 
change is simply given by In (kt + c). The effects of increasing 
k are similar to those for the other two growth laws — increasing 
(AW74) max and increasing the final rate of weight loss. The 
ratio C m /C 0 varies from 0.33 to 0.43 as k varies from 0.050 to 
0.001 and is thus generally higher than the corresponding ratio 
for any parabolic case. 

Figure 5(b) illustrates the effects of increasing m, with k fixed 
at 0.01 . As in the power grow th law case, increasing m reduces 
the final weight loss rate but now also produces a clear reduction 
in the maximum weight gain. Furthermore, the times to reach 
maximum or zero weight gain increase, with the ratio 


CJCq ranging from 0.40 (m = 0.5) to 0.33 (m = 4.0), again 
generally higher than the ratio (0.30) for any parabolic case. 

It is generally true, then, that the power growth law and 
logarithmic growth law cases are more versatile in fitting a 
curve because of the additional adjustable parameters. How- 
ever, this versatility implies a more complex process that 
becomes more difficult to confirm on the basis of established 
physical processes and well-known kinetics. 

Spalling case . — Figure 6(a) compares the Monte Carlo, 
discrete segment, bimodal spalling distribution case against the 
uniform thickness spalling geometry. The same k (0.01), Q 0 
(0.002), and At ( 1 hr) are used. Three bimodal models are shown 
for varying numbers of segments, each model with an /? spa]1 
of 1 .0 (spalling to bare metal). At 10 segments one can see the 
wide diversions caused by sporadic spallation events, at 100 
segments only noisy variations from monotonic behavior are 
seen, and at 1000 segments essentially a smooth curve is 
obtained with a very slight noise factor. The maximum weight 
gains are 0.62 mg/cm 2 for the bimodal case and 0.66 mg/cm 2 
for the uniform spall model. These bimodal curves produce a 
final slope that is 40 percent larger than that for the uniform 
spalling case. However, the two spalling cases appear equiva- 
lent for the initial portion of the curve before reaching zero 
weight change. 

Figure 6(b) shows the effect of varying fl spaI} from 0. 1 to 1 .0 
for the same bimodal case with 1000 segments. At /? spa]1 ^ 0. 1 
the bimodal case (line) is equivalent to uniform spalling (sym- 
bols) for this set of COSP conditions. Therefore, it is usually 
most practical to either use the uniform spall model case or the 
extreme bimodal case with # spa u = 1.0. 




Figure 5. — Effect of logarithmic growth kinetics on baseline COSP model curves for uniform spalling of AI2O3 scale. Spall 
constant Qq, 0.002 cm 2 /mg; secondary growth constant c, 1; cycle duration At, 1 hr. (a) Growth rate k varied for growth 
exponent m of 1 . (b) m varied for k of 0.01 . 
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Figure 6. — Comparison of bimodal (Monte Carlo) spall model curves with uniform spalling case for baseline COSP model. 
Parabolic rate constant k p , 0.01 mg 2 /cm 4 hr; spall constant Qq, 0.002 cm 2 /mg; and cycle duration Af, 1 hr. (a) Serrated 
behavior for 10-segment surface and smooth behavior for >1 00-segment surface at depth fraction P S pall of 1 
(b) Convergence to uniform spalling behavior for 1000-segment surface at fl S pall - 0-1 • 


As originally stated in reference 7, the discrete bimodal 
spalling model is generally not needed for adequate modeling of 
a cyclic curve unless the spalling depth extends to bare metal 
(R au= 1.0). Therefore, it is most compelling touse the bimodal 
spalling case when it is well established that the spalling 
geometry is definitely interfacial in nature, as with poorly 
adherent scales, and not described by fracture at some level 
within the scale. 

In summary, it has been well established that the character- 
istic parameters describing the model curves (C m , C Q , C m /C 0 , 
(AW7A) max , and Final slope) all vary in a regular mathematical 
fashion with the input parameters. These relationships have 
been previously analyzed by regression and presented in detail 
(ref. 7) and were recently summarized (ref. 8). 

Basic Examples 

Nickel Aluminides Plus Zirconium 

Figure 7 shows experimental data obtained for the 1 100 and 
1200 °C cyclic oxidation of NiAl doped with 0.1 at.% zirco- 
nium (Zr) by Barrett (refs. 6 and 16). Also showm are COSP 
model curves obtained by iterative model selections for k p and 
Qq, assuming the parabolic growth law and the uniform spalling 
case. At 1 100 and 1 200 °C the corresponding values of k were 
0.0025 and 0.0170 mg 2 /cm 4 hr, respectively, and the values 
of Qq were extremely low, 0.00002 and 0.000095 cm 2 /mg, 
respectively. These two very well behaved, long-time curves 
lend themselves to excellent fits by the COSP model. 

The corresponding k value was determined from 100-hr 
isothermal tests of the same alloy at 1200 °C to be 0.0276 mg 2 / 


cm 4 hr, or nearly 50 percent higher than those predicted by the 
model. In isothermal tests of similar NiAl + Zr alloys the 1100 
and 1200 °C rate constants were 0.00132 and 0.00553 mg"/ 
cm 4 hr, respectively, or more than 50 percent less than those 
predicted by the model (ref. 17). For low oxidation rates this 
degree of variability is not especially uncommon. Although the 
agreement with independent isothermal tests is within an order 
of magnitude, it is not necessarily precise. The physical process 
(parabolic growth) may somehow be indirectly modified by the 
very action of cycling, but it is not clear how this w T ould happen. 



Figure 7. — COSP model (curves) for fitted parabolic growth 
and uniform spall case with 1100 and 1200 °C cyclic 
oxidation data (symbols) for NiAl + 0.1 Zr using 1 -hr 
cycles (refs. 6 and 16). 
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We suggest that some nonideality in the spalling sequence 
(unpredictable variations in Q 0 ) may produce a skewed fit and 
give the appearance of a different growth constant. Although it 
is straightforward to independently confirm the parabolic growth 
kinetics in isothermal thermogravi metric tests, it is much more 
difficult to establish the true extent and geometry of the spalling 
process over a large number of cycles (ref. 2). Thus, some 
caution is appropriate regarding the degree of comparability 
between cyclic modeling and isothermal tests. 

Single-Crystal Superalloys 

SamplesofPWA 1480 were prepared and hydrogen annealed 
(desulfurized) to various sulfur levels in order to produce dif- 
ferent levels of scale adhesion in 1 100 °C cyclic tests (ref. 18). 
Figure 8 shows that at low sulfur (S) levels (<0.2 ppmw) 
excellent scale adhesion was produced and minimal spalling 
occurred over the 1000 1 -Hr cycle tests (sample 20-^6, H 2 
annealed for 100 hr at 1200 °C to 0.14 ppmw S). Here it is 
immediately established that even the best fit to a parabolic 
model, with k = 1 . 2x 1 0” 3 mg 2 /cm 4 hr and C ffJ /C 0 =0.31, cannot 
adequately follow the shape of the actual curv e. However, a 
power growth law model, with m = 4 and k = 0.00075, produces 
CJC 0 = 0.18 and fits the data exceptionally well. This model 
calculates the fraction of scale spalled to be 0.075 percent and 
the retained scale weight to be 1.49 mg/cm 2 after 1000 cycles. 
The improved fit of the power law model versus the parabolic 
growth model is due to the initial period (20 hr) of strong 
transient oxidation. 

Although it is generally assumed that the protective behavior 
of single-crystal superalloys is derived from parabolic growth 
of a-AI 2 0 3 , it is also known that the initial transient growth of 
NiAl 2 0 4 is relatively large (ref 19), giving the appearance of 
nonparabolic growth. Independent 1 100 °C isothermal tests 
have identified deviations from parabolic growth for approxi- 





Figure 8. — COSP model power law fit to 1 1 00 °C cyclic 
oxidation behavior of PWA 1480 using 1-hr cycles; H 2 
annealed to 0.14 ppmw S (ref. 17). 


mately the first 16 hr of oxidation, producing essentially a 
0.1 1 -mg/cm 2 transient oxidation “y intercept” on a f 1/2 para- 
bolic plot (ref. 20). The parabolic kinetics were described by a 
k of 6.5x1 Cri 4 mg 2 /cm 4 hr in contrast to 1 .2xl0 -3 mg 2 /cm 4 hr in 
tne parabolic fit here. Both parabolic and power growth law 
COSP models yield spall constants Q 0 very close to 
0.0005 cm 2 /mg, which is a relatively small factor owing to the 
strong scale adhesion observed in the low-sulfur samples. 

At a slightly higher sulfur level (sample 20-5, H 2 annealed 
for 8 hr at 1 200 °C to 0.8 ppmw S), more spalling took place and 
negative weight changes were sustained. Figure 9 shows the 
experimental data along with three model curves. The power 
growth law model used the same growth parameters as deter- 
mined in figure 8 for the more adherent scale. The parabolic 
model was best fit by a k of 5xl0 -3 mg 2 /cm 4 hr. The spalling 
constants were 0.010 to 0.015 cm 2 /mg, or >20 times the 
constants determined for the low-sulfur sample in figure 8. 
Some deviation from a good fit Is apparent for the longer test 
times. Here the bimodal case, with /? spaJ1 = 1 .0 for bare-metal 
spalling, shows perhaps the best overall agreement as A1 2 0 3 
interfacial spallation becomes a more prominent feature of the 
cyclic oxidation process (now with = 1 .2 percent, in contrast 

to 0.075 percent for sample 20-6 with 0. 14 ppm S at the end of 
1000 cycles). The calculated weight of the scale retained after 
1000 cycles is 0.93 mg/cm 2 . 

Figure 10 illustrates even higher sulfur samples 
(sample 50-3, H 2 annealed for 20 hr at 1 200 °C to 1 .5 ppmw S; 
and sample 20-0, with 6.7 ppmw S, not annealed). Using the 
same power law growth parameters for AI 2 0 3 scales, the values 
of £>o were now determined to be at the higher levels of 0.05 and 
0.5 cm 2 /mg, respectively. However, these fits were not adequate, 
especially at the longertimes. X-ray diffraction identified nickel 
tantalate(NiTa 2 0 6 ,oran equivalent rutile structure) as the major 
scale phase owing to massive spallation, aluminum depletion, 
and a shift in oxide type from predominantly a- A1 2 0 3 scales to 



Number of 1 -hr cycles, C 

Figure 9. — COSP power law growth and bimodal spalling 
model fit to 1100 °C cyclic oxidation behavior of PWA 1480 
using 1-hr cycles; H 2 annealed to 0.8 ppmw S (ref. 17). 
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Figure 10. — COSP model fit to 1100 °C cyclic oxidation 
behavior of PWA 1480 using 1-hr cycles; H 2 annealed to 
1 .5 ppmw S or 6.7 ppmw S for no annealing (ref. 1 7). 
(Curves denote model predictions; symbols denote 
experimental data.) 



Figure 11 . — Metal consumption rates during cyclic oxida- 
tion of PWA 1480 as function of alloy bulk sulfur content. 


heavier Ni,Ta scales. Thus, for sample 50-3, with 1 .5 ppmw S, 
an estimated stoichiometric constant (from twice as much rutile 
as spinel plus alumina) of 4.3 1 52 is derived, and an excellent fit 
is obtained by a parabolic COSP model. Parabolic fits employ 
heavier nonprotective oxides owing to excessive initial interfa- 
cial spalling. Similarly, for sample 20-0, at 6.7 ppmw S, the 
diffraction data suggest a major rutile phase with lesser amounts 
of a mixed spinel, yielding a stoichiometric constant of 4.8238. 
(By comparison, S c = 5.38 13 for pure NiTa 2 0 6 ). However, the 
multiple inflections in these data indicate numerous mecha- 
nism changes that prevent a fit by just one model. The calcu- 
lated fraction of oxide spalled and the final retained scale weight 



Figure 1 2. — Rough model fit, using mixed oxide stoichiom- 
etry and relatively high growth constant, to erratic 
behavior of Ti-24Ai-1 1 Nb (at.%) in 982 °C cyclic oxida- 
tion with 1-hr cycles. (Curve denotes model prediction; 
symbols denote experimental data.) 

were 0.97 percent and 3.83 mg/cm 2 , respectively, for sample 
50-3 and 2.63 percent and 2.56 mg/cm 2 for sample 20-0. 

The effect of sulfur content on weight change behavior can 
be connected to various model parameters. Initially, Q 0 would 
seem to be the logical figure of merit. However, because of 
changes in the oxide phase and the growth models above 

1 ppmw sulfur, a simple relationship was not found. Instead, 
the final slope of the model curves was examined, as this is 
equivalent to the rate of metal consumed and represents the 
ultimate measure of material wastage. The results plotted in 
figure 1 1 indicate a relatively smooth curve. An increase of two 
orders of magnitude in metal consumption rate occurs between 
0 and 10 ppmw sulfur, with most of the effect between 0 and 

2 ppmw. 

Rather than depend only on the weight change data for these 
mechanistic transitions, it would also be advisable to use models 
of solute depletion and diffusion modeling to define a critical 
flux required to sustain protective scale formation. Such a 
computer model has been developed (COSIM) and coupled 
with COSP to account for these transitions to less protective 
behavior (ref. 21). 

Titanium Aluminides 

Another nonprotective system, with a considerable amount 
of fast-growing titanium dioxide (Ti0 2 ) making up the scale, is 
now considered. High-temperature structural titanium 
aluminides have been engineered with excellent specific strength 
at intermediate temperatures, but they suffer from rapid oxida- 
tion kinetics at high temperatures. Figure 1 2 (ref. 22) shows the 
982 °C cyclic oxidation data for one such alloy, o^-Ti^Al + Nb 
(24Al-llNb in at.%). Considerable weight loss has been 
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sustained after only 100 1 -hr cycles of testing. The data are 
rather erratic, so that no single model can provide a satisfactory 
fit. However, the average behavior has been approximated by 
a model using a 50-50 molar mix of TiC^-A^C^, a k of 

0. 234. mg 2 /cm 4 hr (obtained from independent isothermal tests 
of similar alloys, ref. 23), and a Q 0 of 0.005. By assuming an 
oxide stoichiometry of either pure Ti0 2 or pure A1 2 0 3 , essen- 
tially the same oxidation curve was maintained, with corre- 
sponding Qq of either 0.0035 or 0.0065, respectively. The point 
is that information can be obtained for nonprotecdve complex 
scales with some consistency from independent tests. How- 
ever, the danger always remains that the counteradjustment of 
two parameters may produce a reasonable fit that may not be 
firmly justified by these tests. 

Model Restrictions 

Certain combinations of input variables may preclude a 
physically meaningful or mathematically logical calculation of 
the oxidation curve: 

1. The spall constant input Q 0 shall not be so high as to 
produce a spall fraction F s greater than unity. The 

fraction spalled F s of the existing scale weight W' is given 
by F S = Q 0 ( W' ) a . Therefore, as W' increases with time (as 
governed by k p ), F s also increases. It is then quite possible 
that before the limiting oxide weight is achieved according 
to COSP, the quantity Q 0 ( W' ) a will exceed unity, a physi- 
cal impossibility. At this juncture the calculations are halted 
and a warning message is triggered, recommending that a 
lower value of Qq be tried. 

2. In Monte Carlo cases the input value of # spa u shall not 
be so low as to cause the spall fraction F s to be less than 
that dictated by £> 0 . In the present version of COSP the 
bimodal case of nonuniform, random (Monte Carlo) spalling 
of individual segments is defined by the average spall 
fraction F avg as given by F = PF spall , where P and F spall 
are the probability and specified depth fraction of spalling. 
If # spa ji is set too low, then even with P = 1 .0 (the maximum 
probability), F determined from F avg = P^ spa n may still 
be less than that required by F avg = Q 0 ( W' ) a . If this 
condition occurs, calculations are halted and a warning 
message will be triggered, recommending that a larger 
value of pa]l be tried, (Conversely, a smaller value of Q 0 
may also remedy this inconsistency.) 

3. Multiplicity/ambiguity of input parameters may occur 
in fitting experimental curves. This restriction is a general 
warning that various combinations of input parameters 
may, especially in extreme cases, provide seemingly equally 
acceptable fits to a desired cyclic oxidation curve, but with 
unacceptable physical basis in reality. For example, a strongly 
negative slope in weight change may be provided by the 
product of many combinations of k p and Q 0 . However, only 


a restricted range of k is realistic for a given oxide type at 
a given temperature. Typical values may be provided from 
the literature or from independent tests. It is generally 
advisable to perform sensitivity tests on the effect of input 
parameters on the various descriptive parameters. The spe- 
cific relationship of the descriptive parameters to the vari- 
ous combinations ofinput parameters indicates which factors 
will produce these equivalent multiplicity effects (ref. 7). 

Concluding Remarks 

The preceding text has provided a tour through the major 
operating features of COSP for Windows (WinCOSP), a 
Windows-based computer program developed to reproduce 
and model the cyclic oxidation weight change behavior of 
metals and alloys. Model curves have been presented to show 
the major trends in weight change, scale growth, spallation, and 
material consumption with variations of input parameters for a 
number of growth and spalling cases. For a single well-behaved 
cyclic curve it is especially useful to be able to extract param- 
eters such as the rate constant, the scale thickness, the fraction 
of oxide spalled, and the total amount of metal consumed. 

Examples have been presented to illustrate classic mono- 
tonic behavior and good agreement, such as for NiAl + Zr. Also, 
nonparabolic (logarithmic) cases with large transients and 
deviation from a single spalling model have been proposed for 
the superalloy PWA 1 480. A major discontinuity was observed 
for the higher sulfur contents, with the eventual change from 
primarily A1 2 0 3 to Ni,Ta-rich oxides. Here the use of a higher 
stoichiometric constant markedly improved the model fit to the 
data. Despite these specific idiosyncrasies, the overall model- 
ing effort was informative in determining metal consumption 
rates as a function of sulfur content. 

A number of concerns exist regarding overselling the agree- 
ment of any specific model with a set of experimental data. 
There is never any guarantee that good agreement with one 
model automatically precludes equally good agreement with 
some other model. Furthermore, agreement with 1 000-hr data 
does not necessarily guarantee unqualified predictive ability 
for 10 000-hr data. It is strongly recommended that model 
parameters be critically evaluated by independent measures of 
oxide type, scale growth, and spalling modes. The models 
should be presented as reasonable fits to the data for the purpose 
of showing major trends and consistent near-term predictions, 
but should not be claimed as categorical proof of detailed 
mechanisms for indefinite exposure times. 

COSP for Windows can be obtained for individual use by 
contacting the authors, James.L.Smialek@grc.nasa.gov or 
Judith.V.Auping@grc.nasa.gov. 

Glenn Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, November 2001 
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